Transient dynamics of open quantum systems 
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We present a renormalization group (RG) method which allows for an analytical study of the 
transient dynamics of open quantum systems on all time scales. Whereas oscillation frequencies 
and decay rates of exponential time evolution follow from the fixed point positions, the long-time 
behavior of pre-exponential functions is related to the scaling behavior around the fixed points. We 
show that certain terms of the RG flow are only cut off by inverse time, which leads to a difference 
between infrared and ultraviolet scaling. An evaluation for the Ohmic spin boson model at weak 
damping reveals significant deviations from previous predictions in the long-time regime. We propose 
that weak-coupling problems for stationary quantities can in principle turn into strong-coupling ones 
for the determination of the long-time behavior. 

PACS numbers: 05.10.Cc, 05.30.-d, 05.30.Jp, 73.23.-b 



The time dynamics of a small strongly interacting 
quantum system coupled to noninteracting large reser- 
voirs is a fundamental issue in nonequilibrium statistical 
mechanics. The prototype is a two-level system coupled 
to an environment ji* 2 - which is of particular interest in 
quantum information processing.— A typical setup is the 
one of transient dynamics: The system and environment 
are decoupled for times t < and the coupling is switched 
on suddenly at t = 0. The time evolution of the reduced 
density matrix p t of the local system will then be char- 
acterized for t > by a series of terms, each of which 
will generically be of the form of an exponential together 
with a pre-exponential function p t — J2n -^t ex P {~^ z nt)- 
Here, z n — £1„ — iT n consists of an oscillation frequency 
fi n and a decay rate T n , where one of the scales z n will 
be zero characterizing the stationary state p s t for t — > oo. 
Besides the calculation of p s t and z n , the main challenge 
lies in the analysis of the pre-exponential functions F™ 
on all time scales. Although interesting field-theoretical 4 - 
and numerical techniques 5 - have been developed to study 
the time dynamics, the precise form of pre-exponential 
functions has not been addressed so far. Promising tools 
for this purpose are perturbative renormalization group 
(RG) methods for nonequilibrium problems, like the flow 
equation method,— real-time RG (RTRG)^~— and func- 
tional RG 12 techniques, or combinations of the latter 
two^ The RTRG method allows for an analytical study 
on all time scales, provided that the RG flow stays in 
the weak-coupling regime. The time dynamics is related 
to the density matrix p(E) in Laplace space, where the 
exponential scales z n are the singularities of p(E) in the 
complex plane and the pre-exponential functions can be 
determined from branch cut integrals starting at these 
singularities. In Ref. [ll| a RG approach has been pro- 
posed by using the Laplace variable E itself as flow pa- 
rameter (called E-RTRG in the following) , where the sin- 
gularities z n are given by the fixed points of the RG flow 
and the long-time behavior of pre-exponential functions 
can be related to the scaling behavior around the fixed 
points. 



In this Rapid Communication we will combine E- 
RTRG with a new parametrization of the effective Liou- 
villian in terms of slowly varying logarithmic functions 
and provide a discussion of the generic time evolution on 
all time scales. The main result is the insight that, for 
the determination of pre-exponential functions, certain 
terms of the RG flow are only cut off by the energy scale 
of inverse time 1/t. This is in contrast to stationary 
quantities, where it has been propose d 14 ' 15 and micro- 
scopically show n 16 i 17 that all terms of the RG flow are 
cut off by decay rates. As a consequence, we find that 
the long-time behavior is generically quite different from 
that discussed in Refs.l9land[l0lat intermediate and short 
times. To show this explicitly we will apply our method 
to the Ohmic spin boson model at weak damping, which 
turns out to be a weak-coupling problem even close to the 
fixed points. For the diagonal components of the density 
matrix, we find that the power-law exponent for the scal- 
ing behavior of the pre-exponential function agrees with 
perturbation theory, in contrast with that predicted by 
the noninteracting blip approximation (NIB A). 1,2 For the 
nondiagonal elements we find a rather complex scaling 
behavior which differs from that of perturbation theory. 
We expect similar deviations to occur for other models 
of open quantum systems as well. In particular, for cer- 
tain problems, e.g., the antiferromagnetic nonequilibrium 
Kondo model at large bias voltage, it may even turn out 
that the renormalized vertices are small for the calcula- 
tion of stationary quantities, but become large close to 
the fixed points z n , i.e., a weak-coupling problem for sta- 
tionary quantities can turn into a strong-coupling one 
for the study of the long-time behavior. This opens up 
another class of interesting problems for the future. 

Generic discussion. — For a discrete quantum sys- 
tem coupled to noninteracting reservoirs at time t = 0, 
it can be shown^ that the time evolution of the re- 
duced density matrix p t of the local system follows from 
ipt = J Q dt' ' Lt-t' Pt' i where L t is an effective Liouvillian 
superoperator acting on local operators. Defining the 
Laplace transform via A(E) = J °° dte lEt A t , one obtains 
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the formal solution p(E) — iH(E)p t= o with the propaga- 
tor H(E) = E _\^ E \ j and the time evolution follows from 
the inverse Laplace transform 



Z7T 



oo+iCT 



dEe- tEt n{E) p t=Q 



(1) 



Reference |a shows how to calculate the Liouvillian L(E) 
from a diagrammatic expansion in some appropriately 
defined dimensionless system-bath coupling a. The prop- 
agator II(-E) is an analytic function in the upper half of 
the complex plane with poles at z k (k = st, 0, ±1, ±2, . . . ) 
in the lower half of the complex plane, k = st de- 
notes the zero pole z^ t = 0, which defines the station- 
ary state. At zero temperature^ additional nonanalytic 
features arise from branch cuts starting at the singular- 



plane. Using the general expressions ([T]), Q and ((4]), the 
typical time dynamics can be obtained as follows: 

For short times, t <C l/|z„|, one can replace E — > 1/t 
in the logarithmic functions La and Z'. This gives for 
Eq. © the result p t « e^ 1 ^ 1 /^ 1 Z' (l/t)p t=Q . Expand- 
ing the exponential, one finds in leading order that the 
scaling behavior at small times follows from the scaling 
behavior of Z'{l/t) at large energies. This is the poor 
man scaling regime, where the cutoff scales z n are unim- 
portant. As a result, one obtains universal short-time 
behavior, which has been reported, e.g., for spin bosoni^ 
and Kondo models^. 

For intermediate and long times, t > l/|z„|, we con- 
sider each integral separately around the branch cut at 
E = z n — ix, x > 0. It leads to an exponential factor 
multiplied by the pre-exponential function Fp 



ities z n — z v k + Ap n , which are generically given by the ^ J Q °° e xt Sf(z n -ix), where Sf(E) = f(E+0 + )-f(E- 



pole positions shifted by some linear combination Ap n 
of the chemical potentials of the reservoirs. They arise 
since the Liouvillian depends logarithmically on terms 
~ aln( „_ ); generated by ultraviolet divergencies in 
the band width D of the reservoirs. The logarithmic 
divergencies can be systematically resummed by using 
E-RTRG. 11 The RG equations express derivatives of the 
Liouvillian L(E) by a diagrammatic series in terms of 
effective vertices, which is free of divergences and can 
be systematically truncated for weak-coupling problems. 
Solving the RG equations along the path E — z„+iA±0 + 
in the complex plane starting at A = D, one can tune the 
positions of the branch cuts to E = z n — ix, x > 0. 

For a discussion of the generic time evolution, it is very 
helpful to use the form 



0~) denotes the jump of the integrand across the branch 
cut. We start with the contribution from a branching pole 
z n = zt. Here, one can approximately replace x — > 1/t 
in all logarithmic functions and take the average value 
Pk(z k — i/t) and Z'(z k — i/t) across the branch cut. The 
jump over the branch cut is dominated by a delta func- 
tion from the resolvent $ E-\ k (E) ~ 27T(5(x). This gives 

the contribution p^' p from the branching pole z k 
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(5) 



L(E) = L A (E) + E L'(E) , 



(2) 



where La(E) and L'(E) are slowly varying logarithmic 
functions. This form is valid in the universal regime 
\E\ <C D. The nonuniversal regime \E\ > D, which cor- 
responds to the ultrashort time regime t < 1/D, is not of 
interest here since it depends on the microscopic details 
of the high-energy cutoff function. In the supplementary 
material^ it is shown how the E-RTRG method can be 
used to obtain RG equations for La(E) = Z'(E)La{E) 
and the Z- factor superoperator Z'(E) = 1/[1 — L'(E)]. 
With these quantities one can express the time evolution 
([T]) in a form where the slowly varying logarithmic parts 
are explicitly shown 



We obtain a logarithmic scaling of the pre-exponential 
function F^' p , which follows from the scaling behavior of 
Pk(E) and Z'(E) around the fixed point z k . For interme- 
diate times, where one can expand the logarithmic scaling 
perturbatively in aln(|z ra |£) <C 1, one obtains the weak- 
coupling expansion of Ref. H|. However, for long times, 
the correct scaling behavior has to be determined from a 
systematic expansion of the full solution of the RG equa- 
tions around the fixed points. For k = st and t — > oo one 
obtains the stationary distribution p st — |x s t(0)) from 
Eq. ©, with L A (0)|x st (0)} = 032i 

For a branching point z n , it is more convenient to start 
from Eq. (fTJ) and write for the jump across the branch 
cut SU = n + ((5i)n_ with II± = n(z„ - ix ± 0+). Using 
II± = Tl±Z'±, neglecting terms of O(SL) 3 , and approxi- 
mating x — > 1/t in the logarithmic functions, we obtain 
after some straightforward manipulations the contribu- 
tion p"' b from the branching point z n : 
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-iEt 



y^_L / dE— —P k (E)Z'(E)p t - , (4) 



k,k' 

SL(z n 



27r JO 

ix) 



dx e 



(z n -iX — \ k ){z n - ix - X k ,) 



—P^Z' n p t=0 , (6) 



X k {E) are the eigenvalues of where A n = X{ Zn - i/t). For times t » 1/ 



La(E) with projectors P k {E). The poles z v k follow from 
z k = AfcOzjp. The integral is performed by closing the 
integration contour 7 in the lower half of the complex 



the argument x ~ 1/t in the denominators can usually 
be neglected^ In this case and for intermediate times, 
where logarithmic scaling is unimportant, we obtain the 
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power law F^' h ~ l/t 1+r , where the exponent r follows 
from the scaling of the jump of the Liouvillian across the 
branch cut SL(z n — ix) ~ x r . We therefore set up an 
equation for ^§B-(z n — ix) and solve this RG equation for 
x — > 0. For example, we find r = for quantum dots in 
the charge fluctuation regime and r = 1 for the Kondo 
model and the Ohmic spin boson models For long times, 
one has to consider in addition the logarithmic correc- 
tions from the RG flow close to the fixed points. 

Ohmic spin boson model. — We now apply our flexi- 
ble method to the Ohmic spin boson model at zero bias 
where a bosonic reservoir H ICS = J2 q 0J q a q a g 1S coupled 
to a two-level system with tunneling A, described by 
the Hamiltonian H = — — cr^. The coupling is given 
by V — ^o~zJ2 q 9q( a q + a i?) with ohmic spectral den- 
sity J(ui) = K^2 q g 2 5(Lo ~ w g) = 2irau>6(uj)p(uj), where 
p(uj) — D 2 /(D 2 + lj 2 ) is some high-energy cutoff func- 
tion. We first summarize our results and compare them 
to previous works. For weak damping a <C 1, we find 
three nonzero poles at zq = —iT and z± = ±f2 — iT/2, 
with the effective tunneling £1 = A(^) a and T = wail. 
In leading order truncation it turns out that no branch- 
ing poles appear. Therefore, according to the general 
expressions ([5]) and (JB]), our results in the intermediate 
and long time regime t > 1 /Q can be written as 
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(7) 



(8) 
(9) 
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T- Pt=0 , 



T+ Pt=0 ■ 



(11) 

(12) 



(At) 2 V 1 
where a = ±, t± = ^(l i a x ), and 

ft = {[l + a\n(nt)][l-\n(l + a\n(nt))]y 2 . (13) 

In Liouville space we ordered the four possible states 

as ++, , H — , — h, where ± denote the two local spin 

states. For two 2 x 2-matrices A and B, we defined the 
4 x 4-matrix 



A® B 



A n B A 12 B 
A21B A22B 
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The results (|7l)- p2)) can also be obtained from Born"" 
and the self-consistent BornSi approximation, but the 
unrenormalized tunneling D, — > A appears and the pre- 
exponential functions can only be calculated up to 0(a), 



i.e., the logarithmic function f t in Eq. (fT2|) is missing. 25 
This correction can only be obtained from a resumma- 
tion of all leading logarithmic divergences at low ener- 
gies. This poses the question of why a similar logarithmic 
correction is not obtained in Eq. (jlll) . If one compares 
our result to the NIBA approximation, 1,2 which discusses 
only the time dynamics of (o~ z )t (for (cr x ,y)t=o — and 
(fz)t=o = 1) and reads 



; cos(£!£) 



2a 



(nty 



(14) 



one finds, besides the missing exponential part e in 
the second termj 2 ^ a different power law exponent than 
the one predicted by our result (jTTJ) . Below we will show 
that there is a subtle reason why all leading logarithmic 
divergencies cancel out in Eq. (fTTj) . which is due to the 
fact that the scaling of the vertex and the Z-factors is 
completely different around the fixed points compared 
with the scaling at high energies. Our results are fur- 
ther substantiated by bare perturbation theory up to 
0(a 2 )2£, which confirms that there are no logarithmic 
terms ~ a 2 ln(fM) in the time dynamics, consistent with 
Eqs. (jnjhCBI [note that f t = 1 + 0(a 2 ln 2 (fii))]. 

We now sketch the derivation of our results. First we 
note that the full effective Liouvillian is decomposed as 
L = L A (E) + EL'(E) + L s , where 



L s = inaA 





1 



is the part which arises from integrating out the symmet- 
ric part of the reservoir contractions. The other parts can 
be parametrized in the following way: 
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z ' E\ .. Z(T 





-iT 

1 



A 

Z_A 



As a consequence, the full propagator can be written as 
n = nZ'(l + L s /E), with ft = 1/(E - L A ). In leading 
order, T_ < a 2 Z-A can be neglected and the cor- 
rection from L s influences only the stationary state © 
and the pole contribution ©. The Liouvillian L A has 
four eigen valu es A s t = 0, Ao = — iT + and, neglecting T_, 
A± = ±iy/ Z-A, with corresponding projectors given by 

PstZ' = f J J ' 7 
P ° Z ' = ' Z ) ; 



P±Z = 2 [±JZZ Z. ' r 



The jump of the Liouvillian across the branch cuts needed 
for the evaluation of ^ is parametrized as 
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In leading order, we obtain the following RG equations 
for Z a , T+ and Sjy^. 



dZ+ <r— * Z+ dZ- 2Z- 



dE 

dT + 
~~dE 



E — A CT dE y E-\ 



(15) 



(16) 



^|±L = 2ir9(x)ag 2 Z_ for E = z a - ix (17) 
= An9{x)ag 2 Z + for E = z - ix (18) 



(19) 



The RG equations are coupled to the renormalization of 
the vertex function g = Z + Z_g 2 . The initial conditions 
at E = iD are Z a = 1, T + = 8^± = 0, and g 2 = 1. 
The solution of these RG equations is very different for 
high and low energies. For energies \E\ 3> O we find 
Z± ps (^) 2q , ~g ps 1, and T+ ps 0. According to our 
general analysis, this gives rise to the universal shorttimc 
behavior 



Pi 



1 



lpt= 



which agrees with previous predictions.^* 2 - In contrast, 
around the fixed points we find a different scaling. Close 



to zq = — «r, we get Z^ 
2a In 



SI \2 
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i47ra(§) 2 a;6'(a;). Close to z CT = <xf2 - iT/2, we 



obtain Z 4 
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Z_ ps (g ) 2 , 



r+ ps T/2 + ioflln(l - a In 



-i(E-z a ) 
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(#) 4 , A and 



+ ioflln(l - 

<57 + (z CT — ix) ps — i27ra(-^) 2 j:0(j:). Inserting these results 
in the general expressions §5§ and ([5]), one obtains the 
results ([7l fT3l) using straightforward algebra. As we can 
see, the most important point is that the vertex function 
g = Z+Z-g 2 is approximately a constant for high ener- 
gies, whereas, at low energies, g is nearly a constant. If 
one makes a mistake and takes 5 = 1 for all energies, one 
would obtain the scaling 5j^(zq — ix) ~ x l ~ 2a 6{x), which 
produces the incorrect NIBA result F° b ~ l/(Vlt) 2 - 2a . 
On the other hand, if one neglects the renormalization of 
g and takes it as a constant, one obtains the correct scal- 
ing for long times but the scaling of the Z-factors at high 
energies will change with the consequence of an incorrect 



description of the dynamics at short times. Therefore, 
for a proper description of the time dynamics on all time 
scales, it is crucial to take the vertex renormalization into 
account. 

Summary and outlook. We have shown that the long- 
time behavior of open quantum systems involves loga- 
rithmic corrections which are generically quite different 
from those at short and intermediate times. For weak- 
coupling problems we propose a perturbative RG method 
in Liouville space with a complex flow parameter where 
these corrections can be calculated from a systematic ex- 
pansion around the fixed points. We applied the the- 
ory to the Ohmic spin boson model and found indeed 
a different long-time behavior than previously predicted. 
Moreover, since certain terms of the RG flow are only 
cut off by the scale 1 jt of the inverse time but not by de- 
cay rates T: it is not guaranteed for all models that the 
RG flow stays in the weak-coupling regime close to the 
fixed points irrespective the size of T. A prominent ex- 
ample of such behavior is the antifcrromagnetic nonequi- 
librium Kondo model at bias voltage V much larger than 
the Kondo temperature Tk, which has been shown to 
be a weak-coupling problem for the determination of sta- 
tionary quantities*^— However, by using the E-RTRG 
method of Ref. [n], it turns out that the renormalized 
vertices become strong for E close to the fixed points 
z n = nV — ir, although T is much larger than the Kondo 
temperature Tk, whereas they stay small for E close to 
zero (which sets the point to determine stationary quan- 
tities). As a consequence, the long-time behavior can not 
be calculated from a weak-coupling analysis but a strong- 
coupling analysis is needed, which goes beyond the per- 
turbative RG method presented in this work and poses 
new interesting problems for the future. Other interest- 
ing situations arise when one of the poles z^ is equal 
to zero and is a branching pole, i.e., when there is no 
exponential decay but a nontrivial preexponential func- 
tion with logarithmic scaling. This happens typically for 
problems with quantum critical behavior or for reservoirs 
with a nonanalytic spectral density in the limit D — > oo, 
e.g., for multichannel Kondo or sub-Ohmic spin boson 
models. Such systems are of particular interest since 
the long-time behavior is no longer suppressed by the 
exponential decay and the logarithmic scaling behavior 
becomes more visible. 

This work was supported by the DFG via FOR 723. 
We thank V. Meden for valuable discussions. 
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Supplementary Material 

We here present technical details of the general derivation of the RG equations and the application to the ohmic 
spin boson model. We start with a short summary of the E-flow scheme of real-time renormalization group (called 
E-RTRG henceforth). 



I. THE E-RTRG METHOD 



In the supplementary material of Rcf. [11] it has been shown how to derive RG equations for the effective Liouvillian 
L(E) by using the Laplace variable E as flow parameter. For quantum dots in the spin/orbital fluctuation regime, 
which are coupled to fermionic reservoirs with a flat d.o.s., like e.g. the Kondo model, the following RG equations 
have been derived in lowest order 

d 2 1 



dE 2 



L ^ = \ = \ J d0j i J dQ 2 f\Cd 1 )f\Cj 2 )G 12 (E)U(E 12 +Ld 12 )G 2 - 1 (E 12 ) (20) 



d 

-^-G 12 (E) = - 
dE y ' 



- (1 2) = - f dCo 3 f(Q 3 )G 13 (E)U(E 13 + u> 3 )G m (E 13 ) -(1h2) (21) 



Here, the filled double-circle represents the effective vertex G 12(E) at zero frequencies = ut 2 = which is defined 
by all connected diagrams with two free reservoir lines. The index 1 = ijaa with r\ = ± characterizes the reservoir 
field operators a+ aa — a^ acJ ,a^ ac , — a aa , where a is the reservoir and a the spin index. The vertex has already 
been averaged over the Keldysh indices and its logarithmic frequency dependence has been neglected in leading 
order. The black horizontal lines between the vertices represent the full effective propagator of the local system 
IIi... n = n(£'i...„ + uJi... n ) with 11(E) = E _\^ E j ■ The indices of the energy argument refer to the green reservoir 
contractions running over this propagator, where, for each contraction, the index has to be taken from the left vertex 
which is connected to this contraction. This are precisely the same indices which have to be used to determine the 
energy argument for the vertex right to this propagator. We have defined i?i... n = E + fti... n = E + p,i + ... Ji n , and 
wi . n = a>i + . . . <D n , with \±i = r)ifi ai and (Dj = rjiUJi. Here, /i Q is the chemical potential of reservoir a and u is the 
single-particle energy of the reservoir state relative to the chemical potential. The green lines connecting the vertices 
represent the reservoir contractions. Including the Keldysh indices they are given by 'yff, (ui, a/) — 5hi8(Cj + uj i )^ p (<D), 
where 7 P (u>) = p'f(p'ui) = f(u>) — | + 1 = ryacr, V = —rj'a'a', and /(w) = e „/T +1 denotes the Fermi function at 

temperature T. Since the cross at each line denotes the frequency derivative only ^7 P ((D) = /'(D) is needed and 
there is no explicit dependence on the Keldysh indices. This is the reason why only the vertices averaged over the 
Keldysh indices appear in the RG equations. However, for the perturbative determination of the initial condition it 
is needed to take the symmetric part ^- of the reservoir contraction into account. Symmetry factors ^ arising from 
the diagrammatic rules (when two vertices are connected by n equivalent lines) are explicitly quoted in (|20p . Finally, 
to calculate the frequency integrals at finite temperature, the following approximate form of the propagator has been 
proposed 

nGE+w) « —L—Z(E) , (22) 
w + x(E) 

where Z(E) = ttttf\ 1S * ne -^-factor superoperator and x(E) = Z(E)(E — L(E)). At zero temperature, this 

approximation is not needed since /'(D) = — 5(u>) for T = 0. 

Obviously, the frequency integrals arc well-defined in the wide band limit D — > 00, i.e. we have obtained universal 
RG equations. This is the reason why two energy derivatives are needed for the RG of the Liouvillian, which is 
typical for problems with spin/orbital fluctuations. If the same formalism is applied to quantum dots in the charge 
fluctuation regime, where the dot states are coupled via tunneling vertices to reservoirs with a flat d.o.s., a single 
derivative is sufficient for convergence and one obtains in leading order the RG equation 



d_ 

dE 



L(E) = I A I = - J du, 1 f(u, 1 )G 1 (E)U(Ex+uj 1 )G l (E 1 ) (23) 



Here, in leading order, the unrenormalized vertices can be taken. 

For the ohmic spin boson model, where a local 2-level system is coupled linearly to an ohmic bosonic bath, the 
situation is different since the reservoir contractions are dressed by the spectral density of the couplings, which is 
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linear in frequency. This means that the frequency integral in (|23l) becomes logarithmically divergent and a second 
derivative w.r.t. the Laplace variable is needed. To show how the reservoir contractions have to be determined we 
consider a bosonic bath H res — ^2 q w q a q a q with uj q > and a linear coupling of the form V = 1^2 q 9q{0'q-\- where 
7 is a generic local operator (i.e. acts only on the states of the local quantum system) and g q denotes the coupling 
between the local system and mode q of the bath. For an ohmic bath, the couplings are characterized by the spectral 
density 

J(uj) = Tr^2g 2 q S(uj - uj q ) = 2tt a u) 6{u) p{u) (24) 
q 

where p{lo) = is some high-energy cutoff function, which here is chosen as a Lorentzian for convenience, a 

is a dimensionless coupling constant characterizing the damping, which is assumed to be small a <C 1. Using the 
notation a^ui) = X] 9 a ')9^( w — u, q)dq-i each reservoir contraction can be expressed via the following average w.r.t. to 
the canonical distribution of the bosonic bath 

(a v (uj)a„'(uj')) = Sn-n' S(u — lo 1 ) r)n(r)(jj) —J(u) (25) 

where (a m a v ' q i) — 8 ri ^ ri i8 qq i'qn{'quj q ) has been used, with the Bose function n(uj) — eU1 /T_ 1 ■ Following Ref. [8], this 
leads to the following contraction between the reservoir field operators in Liouville space 

7^f, (w,w') = 8i\i 5(u) + U) 1 ) — J(j]uj) r\p' n(p'ui) 

7T 

= 6 lV 5(Q + Q') lap' u n(p'u>) %cj) p(u) (26) 

where the indices 1 = rj and 1 = ~r\ involve only the creation/annihilation index and we have used p(u>) = p(—uj). 
Since the local vertex operator 7 is independent of rj, the contraction (|26l) can be averaged over 77 and 77', which gives 

7 pp '(u;,u/) = ^ 7 ff,'(w,w') = 8(lu+lu')Y'{u) , 7 p '(w) = 2ap' u n(p'w) p(w) (27) 
ii' 

In the wide band limit D — > 00 we omit p(cu). Using n(—uj) = —(1 + n(w)), we can split the Bose function via 
n(uj) = — i + (n(u) + ^) in symmetric and antisymmetric part, and obtain for (|27p 

7 p'H = 7 p'H + 7a ( w ) , 7 p'H = -ap'cv , 7 aM = aw(2n(w) + l) (28) 
This gives for the derivatives of the symmetric part 



7s p 'M - ~ap' , ^7f'H = . (29) 
auj aui 



Most importantly, the first derivative is frequency independent and gives no contribution to the RG diagrams (see 
below). This is the reason why only the vertices averaged over the Keldysh indices appear in the RG. For the special 
case of zero temperature, we get n(uj) = —6(—oj), and the antisymmetric part of the contraction together with its 
derivatives reads 

d d 2 
la M = a\uj\ , — 7 a (w) = asign(w) , -r^la{u) = 2aS{uj) . (30) 
aui duj z 

Using this form of the contractions, the universal RG equations follow from the formalism of Ref. [11] in leading order 

as 

52 T,v\ I X X I f , d 2 j a 



dE 2 ' 



L(E) = ^ ~| = duj^G{E)Tl{E + w)G{E) (31) 



-^G(E) = I X I 1 = : - / duj^G{E)n(E + uj)G{E)U{E + uj)G{E) (32) 

Here, the vertex has no further reservoir indices (only for the case of several reservoirs and different coupling 
operators 7 to the reservoirs, the reservoir index has to be retained). Furthermore, the energy argument is just the 
Laplace variable E since there is no chemical potential in the bosonic reservoir (we allow only for energy exchange 
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with the environment). The frequency integrals are convergent and we note that the derivatives of the symmetric 
part of the contraction do not contribute in both RG diagrams. To see this for the RG diagram of the vertex 
renormalization, one closes the integration contour in the upper half of the complex plane and obtains zero since the 
derivative of the symmetric part is frequency independent and the propagator is an analytic function in the upper 
half. However, the symmetric part will give rise to a perturbative correction in the initial condition, see below. 

Decomposition of the Liouvillian. — Due to logarithmic divergencies in the high-energy cutoff D, the Li- 
ouvillian L(E) will generically depend on various logarithmic terms ~ a\n( E L) z ), which are cut off at high energies 
by D and at low energies by the singularities z n of the propagator H(E). In addition, since the Liouvillian has the 
dimension of an energy, linear terms in the Laplace variable E can occur. In the universal regime E <C D, the 
following decomposition is useful, which explicitly exhibits the slowly varying logarithmic parts 

L(E) = L A (E) + EL'(E) . (33) 

Here, L&{E) is proportional to some physical energy scale A (like e.g. temperature, magnetic held, decay rates, 
chemical potentials, level spacing, etc.) but not the Laplace variable E, and La(E) and L'(E) are slowly varying 
logarithmic functions w.r.t. E. With this form we can write the propagator as 

Hi... n = n( J B 1 ... rl +cD 1 ...„) = — l — — (34) 

&1...TI + Wl...n — M-^l-.-n + w l...nJ 



1 



-ZL.n , (35) 



Ei... n + LO\... n — LA(El... n + U>l... n ) — (El... n + Ui... n )L'(Ei,„ n + U>l... n ) E - 

where 

xt..n = <*>l...n + P<l...n — L^. n , L^.. n — Z[ n L^ n , ^i... n ~ \ ' Z'(E) = — JTJj^ > (^6) 

together with 

^i...n = L£,{Ex... n + wx...n) , L l n =L l {E\... n + Q\.„ n ) , Z[ n = Z'{Ex... n + u>i...n) ■ (37) 
Note that Z(E) = pzg~xTE) an< ^ Z'(E) = l _l^ E ) are not identical since 

d E L(E) = L'(E) + {d E L A (E) + Ed E L'(E)} = L\E) + 0(G 2 ) . (38) 
Below we will show that d E L&{E) + Ed E L'(E) ~ 0(G 2 ), i.e. Z(E) and Z'(E) are the same in leading order up to 

We introduce the following graphical notations for the lines connecting the vertices 

ni...„= • • , z[^ n = m- -• , XL..»ni...„= • — - — • , (39) 

such that the relation Z' x n — £JIIi... n = Xi' rJ^i... n can be written diagrammatically as 

-• - E % • = • • (40) 

To find RG equations for L&(E) and L'(E), we consider the three cases of spin/orbital fluctuations (e.g. Kondo 
model, Eq. (|20|)). charge fluctuations (e.g. quantum dots, Eq. (|23|) ) and energy fluctuations (e.g. spin boson model, 
Eq. (|3"T|)) separately. For charge fluctuations, the RG equation (f23| is already of first order, i.e. the Liouvillian 
L(E) « L/\(E) is in leading order a logarithmic function and we get 

Charge fluctuations: ^-L A (E) = - \ * I , J?=L'(E) = (41) 
oE 9 m oE 

For spin/orbital fluctuations we write the second order differential equations (|20j) formally as 
d 2 1 



dE 2 



Thereby, the .E-derivative of the Z-factors and vertices in the second diagram on the r.h.s. gives rise to terms of 
0(G 3 ) and are added for convenience. We now identify the r.h.s. of this equation with the expression 

d 2 d did d 1 



dE 2 v ' dE dE 1 dE dE 



9 



and find 



■S-L'(E) = - H 9. il (44) 
dE y ' 2 14 *4 v ' 



i^+^ = >-^ = 5£?3 + 0(G3) ■ (45) 

The second equation shows that the difference between 8eL{E) and L'(E) is indeed of 0(G 2 ) as stated above after 
Eq. (|3"8"]). Solving (g5} for d E L A (E) and using flU) and (gDJ) we find in leading order 



d 1 /i 

Spin/orbital fluctuations : g^L A (E) = - 7 X ^ , -^ L '( E ) = \ X ^ ( 46 ) 



9 1 



x 

As required we find that L A (E) is proportional to some physical energy scale appearing in n - Note that, at finite 
temperature, also the frequencies u>x... n appear in x%, n i which gives a term proportional to temperature for L A (E). 
The same procedure can be applied to the RG equation (|3"TT) for the case of energy fluctuations with the result 

Energy fluctuations : J^ L *( E ) = ' Je L ' {E) = i^^i (47) 

The RG equations (|4l"|) . (|46|) and (|47|). together with the vertex renormalization (|2"Tj) and (|32|) are the final differential 
equations to be solved to obtain the Liouvillian L(E) in the form (|33|). We note that the choice of the second term in 
(1421) seems at first sight not unique since it is of higher order. However, it is not important whether the construction 
of the RG equations for L A (E) and L'(E) is unique at a certain truncation order, but the crucial point is that all 
corrections to this construction can be shown to be beyond leading order. 

With the quantities L A (E) and L'(E) the propagator n(£') can be expressed as 

n(J5) = l —— T = IL(E)Z'(E) , 11(E) = 1 (48) 

v ' E - L(E) y ' v > ' v ' E-L A {E) 

with 

L A (E) = Z'(E)L A (E) , Z'(E) = - _ 1 . (49) 

To evaluate the RG equations one needs the frequency dependence of the propagator Tl(E + This is approximated 
by neglecting it in leading order within the logarithmic parts L A and Z' , i.e. we use 

n(g + ^) « — 1 Z\E) . (50) 

E + uj - L A (E) 

Finally, due to the properties TrL(E) = and L(E) C = —L(—E*) (with the c-transform defined by (L c ) SlS2 , s , s , = 
(L S2Sl .s lS , )*) [8], we note the following properties for the quantities L A (E) and Z'(E) 

TtL a (E) = , TtZ'(E) - 1 , L A {Ef = —L A (—E*) , Z'(E) C = Z'(-E*) . (51) 

Time evolution. — Once the quantities L A (E) and Z'(E) are known, the time evolution of the local density matrix 
pt can be discussed in a straightforward way. Choosing any initial state pt=o at t = 0, we get in Laplace space the 
solution p(E) = J Q dte %Et pt = iH(E)p t = , and in time space via inverse Laplace transform 



o 

Pt = ^Z dEe-* Ef p t=0 = - / dEe~ iEt 1 Z'(E)p t=0 , (52) 

2^7-00+^0+ E — L(E) 2tt 7-oo+io+ E - L A (E) 



The last form of (|52")l is very helpful for the evaluation of the energy integral because it explicitly exhibits the slowly 
varying logarithmic functions L A (E) and Z'iE). To identify the singularities of the integrand we use the spectral 
decomposition of the Liouvillian L A (E) in terms of its eigenvalues Xk(E) and corresponding projectors Pk{E) 

L A (E) =^\ k (E)P k (E) (53) 
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Since we deal with a non-hermitian superoperator, we have to distinguish the left and right eigenvectors, which we 
denote in Dirac notation by \x k (E)) and (x k (E)\ 

L A (E)\x k (E)) = X k (E)\x k (E)) , (x k (E)\L A (E) = {x k (E)\X k (E) . (54) 

The eigenvectors fulfill the orthonormalization condition (x k (E)\x k >(E)) = 5 kk i and the projectors are given by 
P k (E) = \x k (E))(x k (E)\ with E k Pk(E) = 1. 

Due to the condition TiL A (E) — 0, we obtain either X k (E) — or Tr \x k (E)} = 0. Therefore, the Liouvillian has 
always an eigenvalue zero, which we characterize by the index k = st since it corresponds to the stationary state (see 
below). The other eigenvalues are numerated by k = 0, ±1, ±2, We get 

Tr\x 8t (E)) = ^Hz 8t (£)) = 1 , (x st (E)\ss) = 1 (55) 

S 

Tr\x k (E)) = ^2(ss\x k (E)) = , for k = 0, ±1, ±2, . . . . (56) 

s 

As a consequence we get P s t(E) — \x s t(E))Tr and the property TrL A (E) — Tr L'(E) — can also be written as 

P st (E)Z'(E) = P st (E) , P st (E)L A (E) = . (57) 

Due to the condition L(E) C = —L(—E*), the eigenvalues and projectors occur always in pairs (except for k = 0, st 
where we define k = —k) with 

X^ k (E) - -X k (-E*) , P_ fe (£) = P k {-E*Y . (58) 
Using the spectral representation, the time evolution can be written as 

k < 

where 7 is an integration contour which encloses the lower half of the complex plane including the real axis. Poles are 
located at E = z k = X k (z k ), where zf t = is a pole at the origin. At zero temperature, which we consider from now 
on, additional nonanalytic features occur from branch cuts since X k (E), P k {E) and Z'(E) depend logarithmically via 
terms ~ ln( E ® z ) generated by the ultraviolet divergencies from the high-energy cutoff D (at finite temperature the 
branch cuts turn into an infinite number of discrete poles separated by 2ttT). From the structure of the perturbation 
theory it can be seen that the singularities z n are associated with poles of the propagators H(Ei... n ), i.e. are located 
at -&L... n = z v k , where = E + £ti... n . Therefore, the singularities z n — z v k — pLi... n are generically given by the 

poles shifted by some linear combination of the chemical potentials of the reservoirs. 

Using the general expressions (|52")) and (|5^]) . one can discuss the qualitative form of the time evolution in different 
time regimes. For short times t <^ l/|z„|, we obtain E ~ 1/i > \z n \, i.e. the cutoff scales z n in the logarithmic terms 
are unimportant and can be neglected. Furthermore, in leading order, we can replace E — > 1/t in the logarithmic 
parts. This means that we cut off the poor man scaling equations for L A (E) and Z'(E) at the scale E ^ 1/t and 
obtain from (|5"2"]) 

Pt = f / dEe~ iEt 1 Z'(l/t) Pt=0 = e^mt Z ' {l/t)pt=Q . (60 ) 

Expanding the exponential one finds in leading order that the logarithmic scaling of Z'(l/t) at high energies determines 
the short time behavior. 

For intermediate and long times t > l/|z„|, we have to study the contributions from the poles and branch cuts in 
detail. All branch cuts are chosen to point into the direction of the negative imaginary axis, i.e. are located at z n — ix 
with x > 0. This can even be enforced numerically by solving the RG equations along the two paths E — z n + iA ± + 
with A real and initially given by A ~ D. Since no singularities are surrounded by the two paths, the RG flow is 
analytic and can be used to determine the jump across the branch cut. This is a particular advantage of the E-RTRG 
method, which uses a complex flow parameter in Laplace space. Our choice for the direction of the branch cuts is 
very convenient since e~' lEt = e - lz nt e - xt [ s exponentially decaying in xt, which allows an analytical evaluation of the 
branch cut integrals for intermediate and long times. We start with the contributions from the branch cuts starting 
at a pole or branching pole at z k , which we evaluate by using the form (|59[) . For the branch cut integral we set 
E = z v k — ix ± + and replace in leading order X k (E) — )• z k and the logarithmic function P k (E)Z' (E) by its average 
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Pk{Zk — i x )Z'(zP — ix) over the branch cut, where A(E) — \(A(E + + ) + A(E — + )). Furthermore, in leading order, 
we can use x — > l/t in the logarithmic functions. This gives the contribution 

#> « e-^l-[_ dx e- (— VoT-—^) h(4-mZ>(4-mP t = ■ (61) 
Using _ ix \ 0+ - _ ix 1 _ + = 2tt5(x), we obtain 

p** « e "<* P fc (z£ - i/t) Z'(4 - i/t)p t=0 , (62) 

i.e., for z v k — ttk—iTk, an exponential time evolution with oscillation Qk and decay rate T^, modulated by a logarithmic 
scaling function. For the special term k = st, where z% t = 0, P s t(E) = \x s t(E))Ti and P s t(E)Z'(E) = P s t(E), we get 
the following contribution to the time evolution 

p?' p « |5rt(-i/t)> ^> Pst = |^(0)} , (63) 

i.e. we see that for t —> oo one always gets the stationary distribution p st but, if z p t is a branching pole, logarithmic 
corrections can occur for the time evolution which do not decay exponentially. We note that for the models discussed 
here, there is no logarithmic term in the diagrammatic series involving the pole z P t . The reason is that the projector 
P s t gives always a regular contribution, provided that the symmetric part j p (uj) of the contraction (|28p is an analytic 
function [8]. Providing there is no accidental pole z p _^ st = 0, the pole at E = is isolated and has no attached branch 
cuts. 

The evaluation of a branch cut starting at a branching point z n which is not a pole is more subtle since both Xk(E) 
and Pk(E)Z' '(E) can be discontinuous and cancellations can occur between the two contributions. Therefore, it is 
more convenient to start from the first expression of (f52"| involving the propagator H(E). Denoting by 6 A = A + — A- 
the jump across the branch and by A — h(A+ + A-) the average value, with A± = A(E ± + ) = A ± hSA, one finds 
for the jump of the propagator expanding in small SL ~ GL 

5TL(E) = n+SLIl- = — l -^5L — + 0(5L 3 ) . (64) 



E-L E-L 



Using AB - AB — \SASB, we get 



^ E whr PkZ ' + = E -W^hZ' + 0(<5L 2 ) (65) 



E-L E-L v ' ^E-\k ^E-\ k 



Inserting this in (|64j) . neglecting 0(5L i ), and approximating E = z n — ix — > z n — i/t in the logarithmic functions Afe, 
Pfc and Z', we get the following result for the branch cut integral 

p n,b „ e - iZn t J- / dx e~ xt P fe " Z' n SL(z n -ix) _. P$ Z' n p t=0 , (66) 



27T kk ,_£ n Jo z n — ix — z n — ix — AJJ, 

where A^ = Xk{z n — i/t), Pj? — Pk(z n — i/t) and Z' n = Z'(z n — i/t). We have omitted the cases k = n or k' — n since 
we consider a branching point and not a branching pole. Since AJ! ~ 0(zT), we can neglect x in the denominators of 
the resolvents for times t ~ 1/x 3> l/\z n — z^ k , \. For special resonant cases, where z n comes close to z k or z k , , one 
can also define time regimes l/\z n \ t <C l/|z„ — fc , |, where cc dominates in the denominators for certain values of 
k or k! . In any case, to evaluate the integral over x, it is necessary to know the jump of the Liouvillian SL(z n — ix), 
for which we will derive RG equations in the following. 

RG equation for SL. — In leading order, the jump SL of the Liouvillian at a branch cut with E = z n — ix ± + is 
generated by some propagator IIi...,, in the perturbative expansion, which is resonant, i.e. the jump of this propagator 
across the branch cut becomes a J-function. To tune the branch cut of each propagator w.r.t. E along the direction 
of the negative imaginary axis, we first close all integration contours over the real frequencies Co in the upper half of 
the complex plane, where the only nonanalytical properties are those of the Fcrmi/Bose- functions on the positive 
imaginary axis. This turns all frequency integrations to ones along the positive imaginary axis J dui — > i J °° dw and 
the sign- functions of the antisymmetric part of the Fermi/Bose functions have to be replaced by their jump on the 
imaginary axis sign(w) — > 2. A particular resolvent containing the eigenvalue Xk{E\... n + i&i...n) wu l then become 
resonant if the condition z n = z k — pi... n is fulfilled. With Ex,,_ n = z n + p,\... n — ix ±0 + — z k — ix ± + , we replace 
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approximately X k — > z k , P k — > Pk{z k — ix) and Z' — > Z'(z k — ix), which gives for the jump of the propagator the 
following (5-function 

SRi n = ( t" r ^ hizl - ix) Z'(zl - ix) 

\-ix + iu! 1 ... n + 0+ -ix + iu>!... n - 0+ J Kk ' v k ' 

= 27r<5(u>i... n — x)Pk(zl — ix) Z'(z% — ix) . (67) 

Since u)\... n > 0, the frequency integrals give only a contribution for x > 0. Diagrammatically, we indicate the jump 
of the propagator by 

The RG equations Q^r{z n ~ix) can be obtained by 
We obtain in leading order 

Charge fluctuations: 

Spin/orbital fluctuations: 



Energy fluctuations: 

We have chosen the number of derivatives by the criterion that the frequency integrals on the r.h.s. are convergent. 

These equations are nontrivial only on the branches, where E = z n — ix and z n — z v k — fii... n . Setting this energy 
argument in the propagators between the vertices, we explicitly obtain together with (|6T|) (note that we consider zero 
temperature) 



2ni 0(x) G 1 {z p k — pi — ix) P k {z p k - ix) Z'(z p k - ix) Gi(z% - ix) (72) 



-2tt6(x) G 12 {z p k - (ii2 - ix) Pk{z p k - ix) Z'(z k - ix) G^\{z v k - ix) (73) 



— 4ira9(x) G(z n — ix) Pk(z n — ix) Z'{z n — ix) G(z n — ix) (74) 

where z n = z k in the case of energy fluctuations, and we have replaced all vertices by their average G across the 
branch cut in case that they are discontinuous. The initial condition for the last two equations is 8L(z n ) = 0. 
Up to the corrections from the logarithmic functions, we obtain 8L(z n — ix) ~ 6{x) for charge fluctuations and 
5L(z n — ix) ~ xO{x) for spin/orbital and energy fluctuations. Therefore, if x can be neglected in the resolvents of 
the integrand of ([66)) , we obtain (up to logarithmic corrections) p1' h ~ 1/t for charge fluctuations and p"' b ~ 1/t 2 for 
spin/orbital and energy fluctuations. 



a similar technique as the RG equations for -^L(E) and -^^L{E). 

*HE) = (69) 

^L(E) = (71) 



Charge fluctuations: 

6L(z n — ix) 

Spin/orbital fluctuations: 

d 



,^5L(z n - ix) 



Energy fluctuations: 

d 



. ? SL(z n - ix) 



II. APPLICATION TO THE OHMIC SPIN BOSON MODEL 

We now apply the formalism to the ohmic spin boson model at zero bias, defined by the Hamiltonian H to t = 
H + H rcs + V with 

H = - — (T x , Hres = ^2,U q a\a q , V = ^a z ^2 9q( a q + a t ) • ( 75 ) 

i i 

According to (|24[) . we use an ohmic spectral density and the reservoir contraction in Liouville space is given at zero 
temperature by (|2"8j) 

l p ( w ) = Is (u) + 7o(w) , Is ( w ) = -ap'u , 7o(w) = a\u)\ . (76) 

To set up the algebra for the Liouvillian L(E) and the vertices G(E), we denote the two spin states of the local 
system by ± and order the four states in Liouville space asH — h, , H — , — h This means that states in Liouville space 
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are vectors with 4 elements, corresponding to operators in usual Hilbert space. Superoperators acting in Liouville 
space are 4 x 4-matrices. To parametrize an arbitrary 4 x 4-matrix we decompose it in four 2 x 2-blocks, each of which 
can be decomposed in the basis of the unity matrix do = la and the three Pauli matrices a% (z = 1,2,3 = x, y, z). For 
an arbitrary 4 x 4-matrix A, we introduce the following elegant tensor notation 

(=0,1,3,3 V 1V 11 7 j=0,l,2,3 \ / j=0,l,2,3 

where Ai is a 2 x 2-matrix for all i = 0, 1, 2, 3. This notation has the advantage that a product of two 4 x 4-matrices 
A and £? can be written as 

AB = ^{A 4 ®a i ){B j ®aj) = ^{A i B j )®(a i aj) = {A- E)®1 2 + i{AAB)®a , (78) 

ij ij 

with A T = (A x , A yi A z ), B? = (B x , B yi B z ) and a T — (a x , a y ,a z ). The inverse of a matrix is given by 

a- 1 = Yl A * x ® ■ ( 79 ) 

8=0,1,2,3 

Using this notation we get from the Hamiltonian the following matrix structure for the bare Liouvillian and the 
bare vertices 

L<°> = [#,•] = A J) ®t_ , (80) 
G (o) = £ G pp,(o) = = '0 0\ j (81) 

G (o) ^ P G^) = = (j J W , (82) 

where p is the Keldysh index, [•, •] denotes the commutator and {-, ■} is the anti-commutator. Instead of la and a Xl 
we use the matrices r± defined by 

t± = i(l±(T x ) , with T^r^/ = 8 aa ,T a , r^ 1 = T ff , a z T a = T- a a z . (83) 

Using (|78H83[) . we get for the bare propagator 

arm • n<0> ' + = ■ n(0> - = jE^(^) .(") 

(J V 7 CT=± V ' 

and for the sequence of two vertices 

G <».n<».( E + „) G <». = I E ^-l-^(°j)« T+ + ^-(°;)« T _ , (85) 
G <»)n<»)( E + „ )G <») = i E ^_ s (;°)« T+ . (86 ) 

The vertex can only occur if the symmetric part 7f (lu) = —ap'uj of some contraction connects this vertex with 
a vertex standing left to it. If no other vertex stands between these two vertices, we see from (|86)l that the 
propagator gives a contribution ~ -=y for large u>, i.e. together with the linear frequency dependence of 7J (w) the 
frequency integration J involves an integrand ~ — for large u>. This gives the following contribution to the effective 
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Liouvillian for E = E + i0 + slightly above the real axis 

£~ JL ~[ = -a J dujuj G&> n(°) (E + u + i0+) G^ 
1 f, auj /0 



s 

L s = 



<T = ± 



-crS + A / 



— a / rfw i ., 

2 ^-f 7 S + w-crA + iO+V 10 

~a (-° E + A)(-mt) J duj 5(E + u - oA) ° Q J ® t+ 



<T = ± 

= £tt a A f J ° ^ <g) r + . (87) 

We note that the frequency integral is not logarithmically divergent and can be directly calculated for D — > oo, whereas 
the sequence (|55|) together with the antisymmetric part of the contraction leads to a logarithmically divergent integral 
which has to be treated by RG. The term L s gives rise to a perturbative and energy independent contribution to the 
Liouvillian. It is not possible that the symmetric contraction 7^ (ui) crosses over more than one propagator, since 
otherwise the integrand will be ~ -=j for large oj and the integration contour can be closed in the upper half and 
gives zero since all propagators and the symmetric contraction are analytic in the upper half. Furthermore, since 
Z*Z,(°) = L S G^ = 0, the part L s can not appear in any diagram involving more than two vertices. Therefore, in the 
universal limit D — > 00, we can split the Liouvillian exactly as 

L(E) = L(E) + L s with L(E) = L A (E) + EL'(E) , L s = in a A f ° q) ® T+ ' ^ 

where L A (E) and L'(E) are logarithmic functions which can be determined from a diagrammatic series involving 
only the vertex G^ averaged over the Keldysh indices. Since these diagrams involve always an even number of 

vertices, we find together with H^(E) — J^II^ ® r a , G' ' = ( ^ *j* J ® <r z and the algebra (|8"3")l . the form 

La(E) = Yl^L^E) ® T a and L'(E) = J2<jL' a (E) ® t ct , where each diagram can only contribute to the matrix 

elements L A (E)22 and L' a (E)22- Since the bare quantities are given by = and L''^ = 0, we can parametrize 
L A (E) and L'(E) in the form 

La(E) = ( ° _ ) ® r + + ( ° _ ) g r_ , (89) 

L'(E)=J2(° Q _J im )*T. . (90) 

This gives for L&(E) — Z f (E)L/x{E) and Z ! (E) = 1 _ L 1 /^ the parametrization 

iA(^) = ( J „ lT +{E) ) ® ^ + ( Z _ ( £ )A _ ir A (S) ) ® T " > ( 91 ) 
Z 'W=T{ \ Z?{E))^ ' W 
with Z a {E) = 1+ . r 1 / t^t and ^(E 1 ) = Z CT (J5)r£(.E). Due to L S L(E) — 0, we can write for the full propagator 

U(E) = ^-1^ = U(E) (l + L s I) with ft(£0 = = fi(E) Z'(iS) , (93) 

with 11(E) = — ^ . In the RG equations, only the part 11(E) of the propagator contributes. 

In contrast to the Liouvillian, any diagram for the effective vertex G(E) will involve an odd number of bare vertices, 

which leads to the general form G(E) = J2i= y z 9i(E) ( q *j* J ® involving the sector of the Pauli matrices cr v , z - 
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However, by inspecting the sum of mirrored diagrams, it turns out that only the sector a z remains. To see this 
consider a diagram of the form 

G (0) n (0V G (0) n (0),-a . . . G(°) n^;; ff = A®{a z r a a z r_ CT )» a z = A ® (r_ CT cr.) , (94) 

where X$ is the set of frequency indices crossing over the i-th resolvent from the left, and H^' a = H^^^E + (DxJ- 
Adding the mirrored diagram 

G (o) n (ow G (o) jjcou ^ . . . n (0),-. G( o) n (ov G (o) = A 3 ((Tz T _ ff ^ T£7) „ a x = A <8> (r CT a z ) , (95) 
we get in total A (g) cr z . Therefore, the effective vertex is given by the parametrization 

G(E) = g(E) (° Q ®a z . (96) 

From (|9T]) one can find the four eigenvalues Xk(E) of La(E) together with the projectors Pk(E). A straightforward 
algebra gives the result (we omit the energy argument E in all expressions) 

A st = , P st Z' = ( J ° ^ ® r+ , (97) 
A = -*r+ , P Z' = ( q ^ + ) ® r+ , (98) 

^£_ + ^/z_A2-r 2 _/4 p. z> = -^= <T ( K za' zZx )® T - ■ (99) 



2y / Z_A 2 -T 2 _/i 

With these quantities, the propagator II(i? + cj) appearing in the RG equations can be expressed as 

U(E + U ) « ^— - = ^ F ^ 1 . P*(£) Z'(£0 , (100) 

E + lu - L A (E) E + uj - \k{E) 

where we have used the approximation (|50p neglecting the frequency dependence in all logarithmic functions in leading 
order. 

With the parametrization (I96H100I) for the vertex and the propagator, we now can easily evaluate the RG equations 
(14"T1) and (|32p with Ii(E + w) — > Yi(E + w) and the form (f5T7|) for the derivatives of the antisymmetric part of the 
contraction at zero temperature. We get 

£=L A (E) = / dJ^)G{E) "- L i E) . - Z'{E)G{E) = -2aG(E)-M^- Z\E) G(E) 
oE J duj 2 E + oj - La(E) E-La(E) 

2a E p _ G ^ Pk(E) Z'(E) G(E) , (101) 



OE J da; 2 E + u-L A (E) E-L a (E) 

Xk(E) 
E-X k {E) 

^L'(E) = J*%i»)GW I *WGW = 2aJ2 ¥ -^G { E ) P ki E)Z>(E ) G(E ) 



E + u-L A (E) ' k E ~ x k(E) 

(102) 



8E G ™ = - / *> G(E) g + J £ Z'(E) G(E) g + ^E) G(E) 



-a 

kk 



Y,G[E)Pk{E)Z'{E)G{E)Pk.{E)Z'{E)G[E) f d W sign(o 



1 



P + w - A fe (£) £ + w - X k '(E) ' 

(103) 



and 



= 2, £ Z'(£)G(£) P,(£) Z'(f;)C(J!) ^L^ 9 , (104) 
— Z'(25) = Z'(E) — Z'{E) = 2a £ G(25) P fe (£) Z'{E) G(E) Z'{E) - ^ . (105) 
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Using (we omit the argument E everywhere) 



Z'GP st Z'G = , Z'GPqZ'G = Z+Z-g 2 ( ° Q J ) ® t_ , (106) 



Z'GP a Z'G = Z + Z_g 2 ( ° J ) ®r+ , (107) 



1 



GP fc Z'GP fe ,Z'G = Z + Z_. 9 3 ^= gAg ( ° J UtWwT+liMW-K , (108) 

/ da;signH pT^ p + J-v =- 2 A^v ln l^; ' (109) 

we obtain the RG equations 

^ = mg £ ,/ A2 _ r , /4 ^A7 ' ap =2a5 ^A7 ' (110) 



9 5 - V- ^ 1 E-X 

C-fc , JA 2 -r 2 _/4 A °" A[r E-\ a 

where we have defined g = Z + Z_g 2 , and A = Z-A denotes the renormalized tunneling. The RG equations simplify 
considerably since we can use A ^> r_, which is fulfilled during the whole RG flow since T_ < a 2 A (see below). 
Using A CT w a A and y^A 2 — r 2 _/4 w A, we get the final equations 

sr+ . . ^ Aq - a ct ar_ . r_ - r + 

— — = sao > — , tttt" = 2a o— — , 113 

a 

^± = aZ + ~ g y ^— , d 4- = 2aZ_~g^— , (114) 

dE ^ E — A(7 ' <9£ y £ - A V ' 

CF 

Sir 1 *- '"§^7 ■ < 115 > 



dg _ 1 E - A 

as = A^ ln ^A; ' 

5ff . 2 / 1 1 2 E — Ac 



together with the initial conditions at E = iD 

T±\e=iD = , Z±\e=iD = 1 , g\E=iD — g\E=iD = 1 • (H7) 

To solve the RG equations we first consider the regime of large energies \E\ » |Ao|, |A±|. In this regime we get 

dv+ /r+ z_a 2 \ ar_ r_-r + , N 

~ 2 <*S -F 1 -*-^ > ^«2cg , (118) 



dE y \ E E 2 J ' dE y E 

^F" 2aZ±5 P ' ~dE 3 5 ^3 > (H9) 

a 

which leads to the result 

Z A 2 Z A 2 

T+ w 2ia ~ , T_ « 4ia 2 ~ w 2a T+ , (120) 

-iE\ 2q 2 Z_A 2 



^«(— ■ 5«l-3«^- • (121) 
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Obviously, we have resummed in this solution all powers of logarithmic divergencies ~ (a ln(D/E)) n . This defines 
the poor man scaling regime, which, in time space, corresponds to the short-time regime. Defining a low energy scale 
by fl = A(Q/D) a , we can write the solution for Z± as 



fi\ 2 f-iE\ 2a , . . /f2\ " . /A 




Z , [ -) ) with <> A^-J A 

such that we obtain a universal function in terms of the effective tunneling fl. Since T_ < a 2 A, we disregard it in 
the following. 

We next consider the regime of intermediate energies, where a In rgz^n *C 1 (k = 0,±). In this regime we can 
solve the RG equations perturbatively in a with the result 

^(Ao-A^ln "^"^ , ~g = l + 0(a) , g = (^j + 0(a) , (123) 

!* + «» ""o 1 . ^ lg) 2 (l + 2 Q ln -^ n - Ao) ) , (124) 

where all integration constants have been fixed by comparison with the solution at high energies. In this solution we 
have resummed all powers of logarithmic divergencies ~ (a ln(_D/f2))" , and have expanded in the other small logarith- 
mic functions. It defines the weak-coupling expansion regime, which, in time space, corresponds to the intermediate 
time regime. This weak-coupling expansion is equivalent to the one developed in Ref. [17]. 

From the perturbative solution at intermediate energies, we can already derive the leading order result for the real 
and imaginary parts of the pole positions, defined by zj~ — Xk(zk) (note that Zk = zf. since there is no chemical potential 
in the present problem), with Ao = — iT + and X a = o^J Z-A. For intermediate energies, we can approximately set 
X k (E) « z k (l + 0(a)) in (Tmf and CCM]), and find directly z = -iT + (z ) ~ 0(a) and z a = a y /Z-(z <T )A = 
ail(l + 0(a)). Using this result in (|123[) and (|124[) to calculate the 0(a) correction, we find 

izo = r + (z ) ~ ia^^(z — z a ) In — — ~ —iaVl ^ a ln(icr) = naCl , (125) 

Z -± = y/Z-M * g (l + a hi =fe= Z0) ) * § (1 + a ln(-»)) = g(l ~ i\ a) , (126) 
i.e. with z_ = — zi, the pole positions are given by 

r 

zq = —iT , z± = ±fl — i— with T = 7raf2 . (127) 

Finally, we consider the regime of small energies, where the Laplace variable E approaches one of the singularities 
Zk, such that a In ^jr^ ~ 0(1). First, we note that the RG equation (|115l) for the vertex function g(E) leads to a very 
weak logarithmic correction ~ a(E — z k ) \n(E — z k ) close to the singularities which can be neglected. Therefore, we 
take the constant value g ~ (A/f2) 2 at intermediate and small energies. In contrast, the vertex function g = Z + Z_g 2 
behaves very differently. It is approximately a constant g « 1 for high and intermediate energies, but has strong 
logarithmic corrections ~ ahi(E — z k ) close to the singularities, which arise from corresponding singularities of the 
Z-factors Z±. We start with the fixed point analysis around Zo, where T+(E) ~ T and Z+(E) s» (0,/A) 2 behave 
smoothly. In contrast, Z-(E) has a logarithmic singularity, which can be determined from the RG equation (|114| 

dZ_ 1 , , 1 , /A\ 2 1 , x 

-OE * 2aZ -*— = 2aZ ~ Z + 9 — Q * 2aZ - U J (128) 
1 ,. . o f A \\ -i{E-zo) 



- z_(E) - const - 2a UJ ln ^r^ • ^ 

Fixing the integration constant by comparing with the solution (|124[) at intermediate energies, we find for E close to 
A\ 2 „ „ „ / \ 2 „ fn\ 2 1 



, r + (^)«r , z + (E)*{-) , i 2aln -^ r o) ■ ( 13 °) 
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Close to the singularity z+, we get Z_(E) m (Q /A) 2 but T + (E) and Z + (E) have strong logarithmic corrections, 
which can be determined from the RG equations (|113[) and (|114[) in the following way 

__ „ aZ + .g— - = aZ + Z_g ^— « oZ + (_ j _ - 

1 w const _ Q f^V l n -HE-**) <§P Z + (£0 « f ?V - (P s > ( 131 ) 

sr+ . _a q -z+ . „ „ 2 n . ~ o ESS . 1 n 



iag ~ — iaZ + Z_g as —iaZ + — 



dE E - z + + y E - z+ \n J E - z+ 1 - a ln " t(g ~ z+) E-z 



u 



r + (E) « const + jfllnll-aln-^-^) ^ r+(£) « J + ztt ln f 1 - a In ^ Z 



n ) TV 7 2 v n 

(132) 

where again the integration constants have been fixed by comparison with the solutions at intermediate energies. 
Using a similar analysis close to z- , we obtain for E close to z a the result 

.9(£)«(§) 2 , r + (£) « I + iril In (l - a In * g) ) , (133) 

tt\ 2 1 „ ,„ v /ft x 2 



With these results we can also evaluate the RG equation ([74]) for the jump 5L(zk — ix) of the Liouvillian, which we 
parametrize as 

SL(E) = -i ha{E) ( I \ \ ® r a . (135) 

Using the algebra (fTD5|) and (|TU71) together with our results ((T2"3"|) , ([123]) . ([T5D]) . (fT53"| and (ITM| at intermediate and 
small energies, we get for x < fi from (|74|) 



with the solution 



Q 1 / A \ ^ 

-i— 57+(^(, — ix) = — A-KaQ(x)Z- (z a — ix)g{z a — i^) 2 ^ ~ — 27ra#(x) ( — j , (136) 

-i— — <57_(zo — It) = — 47ra6'(a;)Z + (zo — ix)g(zo — ix) 2 rs — 47ra#(x) ( — ) , (137) 
OS \ SZ / 

/A\ 2 /A\ 2 
5j + (z a ~ ix) « — 27rio; I — J x8(x) , (57_(zo — ix) « —Ania I — J x9(x) . (138) 

Time evolution. — With the results for the Liouvillian we now can evaluate the time evolution. We start with 
the short time regime fit < 1. Using ([5Dj). ([52"]) and (fT2"2l . we find 

Pt ~ ^ q /j^2c J (g) l2/>t=0 • (139) 

The intermediate and long time regime tit > 1 is based on the pole and branch cut contributions given by (|62[) and 
(|66|) . respectively. Thereby, we have to consider that the full propagator (|93|) involves the correction II(_B)L s -g. Since 
L s ~ 0(a), this leads to a negligible O(a 2 )-correction to the branch cut contributions (|66l) but the pole contribution 
(1621) changes to 

P f' P = Art = (Pst(O)Z'(O) + fl(0)L s ) p t=0 « (P st (0)Z'(0) - ]T P fc (0)Z'(0)L s - ] p i=0 , (140) 

\ fc=o,± Zk ) 

p k t p = e-^ t P k (z k )Z\z k ) [\ + L s ^j Pt=0 for fc = 0,± , (141) 
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where we have used P s tZ'L s = 0, Afc(0) « Zk and the fact that all poles are isolated up to leading order truncation. 
Using the form for L s together with the results (I97M99P for the projectors (where we neglect r_), we find 

P st Z'L s = , P Q Z'L S = inaAZ+ ( ® ° ] ® r+ , P a Z'L s = . (142) 



v 1 , 

Inserting these results in (p^CT|) and (fHTj) . we find together with Z + (0) ps Z+(z ) « ^-(z^) w (fi/A) 2 and (TT2"7) 

/ 1 



st,p 



10\ inaAZ + (0) f 0\\ ( 1 0\ 1 



ft =fttK J ^ I 1 ) j T +P t =° ~ I ;^ T +^=° = 2 



1 

n/A 

V ft/A 



(143) 



„'/■" = e-^Z +{zo) [(l ;)+^(; °)}®r + p««e^ ( _„% ( ^ )2 ) r+ ,„-„ . ,144, 



^.P- fi -te.« a (a^/Z-MA Z_(z a )A 



2 1 <tQ/ A (n/A) 2 



T-p t =o ■ (145) 



Finally, by using the algebra of the projectors PkZ' and the jump 5L according to (|97fl99[) and f|135[) . we can write 
the branch cut contribution as 

1 poo „-xt 

PT = e- lZnt ^zZ dx 7 = ufw = P ° 2 '° 5L ^ ^ P °' 2 '° *=° ' ( 146 ) 

2tt ^ J (z -ix- A£) (z - ix - A£, ) 

2^ ^ 7o * - - AJ) 2 



1 Z" 00 p— at 

Pt' = ^ ^zT, dx i — ^ ^ - ix) PS Z"> Pt=0 . (147) 



For Sit > 1 we can use 

zo -ix - A° « -A° w -a^/i^A , z CT - - Aq « z CT 4- i f £ , (148) 
where Z°_ = Z_ (zq — i /t) and = T + (z a — i/t) can be calculated from f|130[) and (| 133[) as 



<> \ 2 1 r 



, FT w - + ia-Q ln(l + a ln(fit)) . (149) 



v A/ l + 2aln(Qi) ' + 2 
Furthermore, due to the algebra of the projectors PkZ' and the jump 5L, we can use 



P P = -t <5 7 + {Z+f [l\j®r + . (151) 
Inserting (|148H151[) in (fl~46|) and ([147)) , and using the result (|1 38[) for Sj-(zq — ia;) and 6~f+ {z a — ix), we obtain 
P?' b « e^ Zot ("2^3) J™ dxe- xt S^(z -ix) (l °) ® r_ Pt=0 



„o-,a -iz„t ( i \ ( z °+ \ r j„ _-irf jc., ^\ ° 



(-#) / dxe x I q J )®r +Pt=0 



«/t Wo o 



(A^ 2 y V 1 



t+ p t =o , (153) 
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where we have defined the logarithmic scaling function 

^ = ((1 + aln(fti))(l - ln(l + aln(m))) 

Eqs. (|143M145l) and (|152M154l) are the final results for the time evolution in the regime Sit 3> 1. We note that the 
branch cut integrals (|146l) and (|147|) can also be calculated exactly in terms of exponential integrals, extending the 
applicability range to the time regime Clt > 1. 



(154) 



